Dynamical polarization of graphene under strain 



F. M. D. Pellegrino, 1 ' 2 ' 3 ' 4 G. G. N. Angilella, 1 ' 2 ' 3 ' 5 ' 4 and R. Pucci 1 ' 2 

1 Dipartimento di Fisica e Astronomia, Universita di Catania, 
Via S. Sofia, 64, 1-95123 Catania, Italy 
2 CNISM, UdR Catania, 1-95123 Catania, Italy 
3 Scuola Superiore di Catania, Universita di Catania, 
Via S. Nullo, 5/i, 1-95123 Catania, Italy 
4 CNR-IMM, Z.I. VIII Strada 5, 1-95121 Catania, Italy 
5 IN FN, Sez. Catania, 1-95123 Catania, Italy 
(Dated: August 27, 2010) 

We study the dependence of the plasmon dispersion relation of graphene on applied uniaxial 
strain. Besides electron correlation at the RPA level, we also include local field effects specific for the 
honeycomb lattice. As a consequence of the two-band character of the electronic band structure, we 
find two distinct plasmon branches. We recover the square-root behavior of the low-energy branch, 
and find a nonmonotonic dependence of the strain-induced modification of its stiffness, as a function 
of the wavevector orientation with respect to applied strain. 
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I. INTRODUCTION 

Graphene is a two-dimensional single layer of car- 
bon atoms, and can be thought therefore as the build- 
ing block of several sp 2 -bonded carbon allotropes, rang- 
ing from three-dimensional graphite, to one-dimensional 
nanotubes, to zero-dimensional fullerenes. Its recent ex- 
perimental fabrication in the laboratory 1 has triggered an 
enormous outburst of both experimental and theoretical 
research. This is justified by the peculiar electronic and 
structural properties of graphene 2 ^, largely due to its re- 
duced dimensionality, as well as to correlation effects. In 
particular, its linear quasiparticle dispersion relation is 
analogous to that of relativistic massless particles, obey- 
ing Dirac-Weyl equation, thus enabling to study quantum 
relativistic effects in a condensed matter system^. 

Most of the unusual electronic properties of graphene 
are encoded in the electron polarizability, which has been 
studied within the Dirac cone approximation at zero 6 
and finite temperature^ for pristine graphene, as well as 
for doped graphene^! These results have been recently 
extended beyond the Dirac cone approximation 10 . The 
effect of spin-orbit interaction in the electronic collective 
excitations of a graphene layer with or without doping 
has also been considered in Ref. [11] 

Here, we will be concerned on the dynamical polar- 
ization of graphene within the full Brillouin zone of the 
honeycomb lattice. While electron correlations will be 
treated at the RPA level, we will explicitly include lo- 
cal field effects (LFE)^, which are characteristic of the 
lattice structure of graphene. The importance of LFE 
have been shown to be more important in graphene than 
in bulk semiconductors, in connection with the static di- 
electric properties of graphene 13 . By discussing the sin- 
gularities of the polarizability, we will be able to identify 
the collective modes of the correlated electron liquid. We 
will be mainly interested in the plasmon modes, which 
dominate the long wavelength charge density fluctua- 



tions. The role of electron-plasmon interaction in renor- 
malizing the (especially low-energy) quasiparticle disper- 
sion relation has been emphasize d 14 * 15 *, and plasmons in 
graphene are potentially interesting for applications in 
nanophotonics^H 



Specifically, we will be interested in the dependence 
of the plasmon modes on applied uniaxial strain. This 
will enable to investigate the interplay between electronic 
and structural properties of graphene. It has been even 
suggested that nanodevices based on graphene could be 
engineered on the basis of the expected strain-induced 
modifications of the deformed graphene sheet (origami 
electronics/^. Indeed, graphene is also characterized 
by quite remarkable mechanical properties. Despite its 
quasi-two-dimensional character, it displays an excep- 
tional tensile strength and stiffness 18 . In particular, re- 
cent ab initio calculations^^ as well as experiments^ 
have demonstrated that graphene can sustain elastic de- 
formations as large as 20%. The possibility of a strain- 
induced semimetal-to-semiconductor transition , w ith the 
opening of a gap, has been therefore studiecP^3 It 
turns out that this critically depends on the direction 
of applied strain, as is also confirmed by studies of the 
strain effect on the optical conductivity of graphene^* 3 -^. 

The paper is organized as follows. In Sec.|TI|we present 
our model, based on a tight-binding description of the 
graphene electronic band structure. We will then derive 
the electronic polarization at RPA level, and explicitly 
include local field effects. We will then derive and dis- 
cuss the various branches of the plasmon modes along a 
symmetry contour of the first Brillouin zone, both nu- 
merically and analytically, in the limit of small wavevec- 
tors. The effect of applied uniaxial strain will then be 
Summary and concluding remarks 



discussed in Sec. Ill 



will be given in Sec. 
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II. MODEL 

A. Tight-binding approximation 

At the tight-binding level of approximation, the Hamil- 
tonian for the graphene honeycomb lattice can be conve- 
niently written as 



H = tia\R)b(R + 5i) +H.c, 



(1) 



Rj 



where a^(R) is a creation operator on the position R of 
the A sublattice, 6(R + Si) is a destruction operator on 
a nearest neighbor (NN) site R + 5^, belonging to the 
B sublattice, and Si are vectors connecting NN sites on 
different sublattices, Si = a(l, \/3)/2, £ 2 = — \/3)/2, 
^3 = cl(— 1,0), with a = 1.42 A, the equilibrium C-C 
distance in a graphene sheeP. In Eq. (51 t £ = t(S e ), 
i = 1,2,3, is the hopping parameter between two NN 
sites. In the absence of strain they reduce to a single 
constant, ti = to, with to = —2.8 eV (Ref. |3T|). 

The dispersion relation of the valence (A = 1) and 
conduction bands (A = 2) are the solutions E^x of the 
generalized eigenvalue problem 



#kUkA = ^kA^kUkA, 



where 



S k 



/ k 
/k o 

1 9k 
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(2) 

(3a) 
(3b) 



and 7k = Ee=ie ik ' Se , h = ELi to*'*' > 5k 
Yle=i S£e lh ' Sl are the usual (complex) structure factor, 
NN hopping, and overlap functions in momentum space, 
respectively. The hopping parameters ti and overlap 
parameters S£ can be expressed in terms of appropri- 
ate pseudoatomic wave functions, which we here take to 
be normalized Gaussian, with standard deviation aj 28 | 32 l 
One finds 



^kA 



-*kTy^+4G k |/ k | 
2G k 



(4) 



where the minus (plus) sign refers to the valence (conduc- 
tion) band, and F k = k /£ + g£fa and G k = 1 - |# k | 2 . 



In the following, we shall also use the abbreviation £ k A = 
^kA — and denote A = 2 for A = 1, and A = 1 for 
A = 2. Moreover, we also set Sa = Si + S2 + ^3 = 0, and 
S b S3. 

A small, albeit nonzero, value of the NN overlap g k 
has the advantage of endowing valence and conduction 
bands with the observed asymmetry. However, since 
g k « 0.077 k <C 1 (also under strain, within the range 
considered below, in Sec. |III[ ), we can safely retain only 
linear corrections to the band dispersions, E\^\ — =p|/ k | — 

+ £K#k)> an d neglect them altogether in the eigenvec- 
tors u kA . 

Our tight-binding approximation is completed by an 
appropriate choice of the Bloch wavefunctions. As in 



zk-R? 



Ref. [El we shall use ^kA = ^" 1/2 Ej <K r ~ ^)Y^' 
where </>(r) is a Gaussian pseudoatomic orbital, and R* 
are vectors of the A = A, B sublattices. 

We can anticipate, at this stage, that some of the find- 
ings of the present study would not be obtained within 
the cone approximation. In particular, the tight-binding 
approximation allows to include important features of the 
electronic band dispersion, such as a finite bandwidth 
and the occurrence of Van Hove singularities. These fea- 
tures will play an essential role in deriving some of the 
characteristics of the plasmon dispersion, which is the 
main goal of the present work. 



B. Local field effects on the electron polarization 

Within linear response theory, plasmon modes can be 
described as poles of the density-density correlation func- 
tion, i.e. the polarization. The random phase approxi- 
mation (RPA) is then the simplest, infinite order, dia- 
grammatic procedure to include electron correlations in 
the dielectric screening giving rise to the polarization^. 
Besides electron-electron correlations, another source of 
k-space dependence of the dielectric function is provided 
by local field effects (LFE) 34 . This is due to the generally 
atomic consistence of matter and, in the case of solids, 
to the periodicity of the crystalline lattice. An account 
of the LFE on the dielectric function of crystalline solids 
dates back at least to the original paper of Adlefl^ (see 
also Refs. I35|36jh and is generalized below to the case of 
graphene, including both valence and conduction bands. 

We start by considering the polarization, which for a 
noninteracting system at finite temperature T reads 



n°(x,x',^ TO ) = h-Hph)- 1 ^^^VkA(x')^A(k,^)V'kA(x)Vk'v(x)^(k',^„ + iw m )^v(x'), (5) 

iuj n kA k'A ; 

where G®(k,iLj n ) = (iu n — <5ca/^) _1 is the Green's function for the noninteracting system, and hu n = (2n + 1)7tA:bT 
[fouj m = 2m7r/cB^l denote the fermionic [bosonic] Matsubara frequencies at temperature T, with h Planck's constant 
and &b Boltzmann's constant. Fourier transforming into momentum space, and performing the summation over the 



3 



Matsubara frequencies, one finds 
n°(q+G,-q'-G',^ TO ) = ^^(q-q')^ £ T kA)k _ qV (i Wm )(k-qA'| e - i ^+ G )-|kA)(kA|e i ^+ G )- f |k-qA'), (6) 

kAA' 



where 



ikA,k-qA'(^mJ — 7^ ^ ? " V 7 J 



6c-qA' — £kA 



Here, tif(uj) is the Fermi function, A c = 3V3a 2 /2 is the 
area of the Wigner-Seitz cell, q, q' belong to the first 
Brillouin zone (1BZ), G, G / are vectors of the reciprocal 
lattice, and LFE are embedded in the Adler's weights^ 



<k-qA'|e- i(q+G) ' f |kA) 

d 2 xe-^^+ G )- x ^kA(x)^ 



k-qA 



(x) 



(-1) 



A — A' 



o i(0 k _ q -0 k )-KJ-<S 3 



-^|q+G| 2 /4 



,(8) 



where in the last line only the onsite overlap between 
pairs of pseudoatomic orbitals, centered on either sub- 
lattices, has been retained, on account of their local- 
ized character, we have retained only the lowest (ze- 
roth) order contributions in the overlap function #k, and 
e z(9 k _ — J k /|J k |. Using a more compact notation, one 
may also write 

II (q+G,-q'-GVa; m ) = {2n) 2 A~ 1 5{q - q') 

x^p qa (G)Q° i8 (q,ia; ro )p^(G') ) (9) 



a(3 



where 



/3* a* (3 
^kA^kA^k-qA'^k-qA' 



N 

kAA' 

X^kA,k-qA'(^m), 



with u^. x the components of UkA (a = 1, 2), and 
p qa (G) = exp(-zG • S a - a 2 g \q + G| 2 /4) 



(10) 



(11) 



are the LFE weights. The continuum limit is recovered 
when G = G' = 0. 

Many-body correlations are then included within RPA, 
yielding a renormalized polarization 

n(q+G,-q'-G',^ m ) = {2tt) 2 A~ x 5{q - q') 

x ^2p (ia (G)Q a0 (q,iuj m )p* cll3 (G'), (12) 



ap 



where now 
Q(q,ioj m ) = 



g s Q (q,ioj m )[l-g s A- 1 V(q)Q (q,icv m )}- 1 , 

(13) 



where matrix products are being understood and g s = 2 
is a factor for spin degeneracy, and 

V a p(q) = X>qa(G")V (q + G" K8 (G") (14) 



is the renormalized Coulomb potential, Vb(q) = 
e 2 / \2eQ£ r q) 1 now a matrix over band indices. Here, 
e r = (e r i + e r 2)/2 denotes te average relative dielec- 
tric constants of the two media surrounding the graphene 
layer. These are air for suspended graphene (e r \ = e r 2 = 
e r = 1). In the case of a stronger dielectric ssubstrate, 
we expect therefore a softening of the correlation effects 
on the plasmon frequency. It is relevant to note that the 
renormalized potential already includes LFE. 



C. Plasmons 

Plasmons are defined as collective excitations of the 
electron liquid corresponding to poles of the retarded po- 
larization, 



n(q,cj) = II(q, -q,zu;„ 



z0 4 



(15) 



where q G 1BZ. Here and in what follows we shall restrict 
to the case G = G' = 0. Indeed, it is apparent from the 
definition of n(q, uS) that its poles can only arise from 
the vanishing of det[H — V(q)Q°(q, ui)\ in Eq. (13), which 
already contains LFE via the renormalized Coulomb po- 
tential, Eq. ( fl4| ) . We therefore define the dispersion re- 
lation Co^(q) of the £-th plasmon branch as 



ir^q.w/tq)) =0. 



(16) 



This clearly involves vanishing of both real and imagi- 
nary parts of the inverse polarization. It will be useful to 
define the dispersion relation d)^(q) of damped plasmons 
through 



Re [n-^q^q))] = 0. 



(17) 



Correspondingly, the inverse lifetime r x (q, uo) of such 
damped plasmons is proportional to — ImII(q, cj), for 
uj = d^(q). 

Fig. [I] shows our numerical results for the plasmon dis- 
persion relation in doped graphene (/i = 1 eV) at finite 
temperature (T = 3 K) along a symmetry contour in 
the 1BZ, without LFE [G" = in Eq. (|l4j, top panel] 
and including LFE (bottom panel). At small wavevectors 
and low frequencies, one recognizes a square-root plas- 
mon mode uj\ (q) ~ y/q, typical of a 2D systenPH This is 
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FIG. 1: (Color online) Plasmon dispersion relation for doped 
graphene (/x = 1 eV) at finite temperature (T = 3 K), not 
including (top panel) and including (bottom panel) LFE. Re- 
sults are shown along a symmetry contour in the 1BZ, with 
r = (0,0), M = (27r/3a,0), and K = (27r/3a, 27r/3>/3a). 
Frequencies uj are in eV. The shaded background is a con- 
tour plot of — ImII(q, lj) (arbitrary scale), while continuous 
lines are the dispersion relation of damped plasmons, u)^(q), 
Eq. (17), is shown as a dotted line. 



in agreement with earlier studies of the dynamical screen- 
ing effects in graphene at RPA level, employing an ap- 
proximate conic dispersion relation for electrons around 
the Dirac point d^P. Such a resu lt has been confirmed 
also for a tight-binding ban a 10 * 37 * , and is here generalized 
with the inclusion of LFE. The effect of spin-orbit inter- 
action can be neglected, in the case of sufficiently large 
chemical potential 11 , as is here the case. 

The high energy (5 — 20 eV) pseudo-plasmon mode, 
extending throughout the whole 1BZ, is rather associ- 
ated with a logarithmic singularity of the bare polariza- 
tion Q°(q, lj) in Eq. (13), and therefore does not corre- 
spond to a true pole of the polarization. This collective 
mode can be related to an interband transition between 
the Van Hove singularities in the valence and conduction 
bands of graphene, and has been identified with a tt — > 7r* 
transitiorPEU. 

At large wavevectors, specifically along the zone 
boundary between the M and the K (Dirac) points, full 



inclusion of LFE determines the appearance of a sec- 
ond, high-frequency (20 — 25 eV), optical-like plasmon 
mode Co>2(q), weakly dispersing as q — >• 0. Multiple plas- 
mon modes are a generic consequence of the possibil- 
ity of interband transitions, whenever several such bands 
are available. This is e.g. the case of quasi-2D quan- 
tum wells (2DQW), whose energy spectrum is character- 
ized by quantized levels in the direction perpendicular 
to the plane of the well, while electrons can roam freely 
within the planed In this case, collective modes arise 
as zeroes of the determinant of the dielectric function. 
At low temperatures, at most the two lowest subbands 
need to be considered. One usually obtains a longitudi- 
nal 'acoustic' mode associated to intrasubband coupling, 
and a transverse 'optical' mode associated to intersub- 
band coupling 3 -^. Such a situation is here paralleled by 
the case of graphene, the role of the two subbands of 
2DQW being here played by the valence and conduction 
bands, touching at the Dirac points in the neutral ma- 
terial. It should be noticed that the plasmon mode due 
to interband coupling is exactly suppressed when LFE 
are neglected. In 2DQW, the discrete nature of the elec- 
tronic subbands is due to the real-space confinement of 
the electron liquid in the direction perpendicular to the 
plane, i.e. to the quasi-2~D character of the quantum well. 
In graphene, the origin of the two bands ultimately lies in 
the specific lattice structure of this material. Therefore, 
the high-energy, 'optical' plasmon mode disappears in the 
absence of LFE (Fig. [I] top panel), as expected when- 
ever the lattice structure of graphene is neglected. In 
other words, while in the absence of LFE only scattering 
processes with momenta within the 1BZ are considered, 
LFE allow to include all scattering processes with arbi- 
trarily low wavelengths, thereby taking into account the 
discrete nature of the crystalline lattice. Such a struc- 
ture needs not be considered in the case of a 2DQW. 
Our finding of a high-energy 'optical' plasmon branch, as 
generic consequence of the two-band electronic structure 
of graphene, should stimulate further investigation of the 
electronic collective modes in graphen d 37 * 40 *, in view of 
the role of electron-electron correlations in interpreting 
electron spectroscopy for interband transitions^. 



D. Asymptotic behaviors 

In certain limiting regimes, one may derive the asymp- 
totic behavior of the polarization in close form. At low 
energies (hu < \t\) and small wavevectors (q — > 0, i.e. 
q<j g <C 1), LFE can be neglected. The matrix product en- 
tering the definition of the polarization through Eq. (13) 
then reduces to 



5 ^ c -V(q)Q°(q, w ) = g.A- 1 V (q)52QP a f l (< l ,u>) 



5 



where Vq = g s (87r/3\^3)(ao/a) Ry, ao being Bohr's ra- 
dius, and 5t(€) = —dn^(e)/de (5(e), as T — > 0. In the 
latter limit, the S- function effectively restricts the inte- 
gration over wavevectors along the Fermi line. Whenever 
the cone approximation holds (i.e., for sufficiently low 
chemical potential and strain; see Sec. Ill), this can be 
taken as the constant-enerev ellipse in Eq~(17) of Ref.EU 
The k-integration in Eq. (18) can then be performed an- 



alytically, and the retarded polarization, Eq. (15), then 
reads 



n(q,cj 



_1 u; 2 g 2 a 2 



where uo^ 



huji 



h 2 ^ 

i0 + , and 



1/2 



|V q £, 



q2 



/a|, 



(20) 



with g s p(jJL) the density of states (DOS) at the Fermi 
level. To leading order in ga, from Eq. (19) one thus 
obtains 



(21) 



for the acoustic-like plasmon dispersion relation. One 
thus recovers the square-root behavior of the plasmon 
dispersion relation, as is typical in 2D electron systems^. 
Moreover, one recovers the dependence of the coefficient 
uj\ ~ n 1 / 4 on the carrier density n, rather than ~ n 1 / 2 , 
as is the cas e for a parabolic dispersion relation of the 
quasiparticles 9 42 . The acoustic-like plasmon mode may 
be related to the Drude weight 43 , thus enabling the ob- 
servation of strain effects from optical measurements^. 
In the case of graphene on a dielectric substrate (e r > 1), 
one has a reduction of uj\ , thus a softening of the plasmon 
mode. From Eq. (19) one may also read off the imagi- 



nary part of the retarded polarization, which close to the 
'acoustic' plasmon mode [u ~ ^i(q)] reads 



Imn(q,w + ) 



'-g^V'^^qa^H^ 



-cJi(q)). 
(22) 

We now turn to the asymptotic behavior of the sec- 
ond branch of the plasmonic spectrum, ^(q). We have 
already established that it displays an optical-like char- 
acter, with Co>2(q) — » ^2(0), as q — » 0. Here, ^(0) is 
greater than the distance between the top of the conduc- 
tion band and the bottom of the valence band. At small 
wavevectors, it is useful to consider the expansions of the 
relevant terms in Eq. (13), which to leading order in qi 
(i = x, y) read 

Qaa(<1,u) « QAA(0,w) + ^2qiyij(w)qj, (23a) 

Qab(<1,u) « -Qaa(0,uj) + ^2q i z ij (uj)q j , (23b) 

ij 

where yij(oj), Zij(uj) are real valued functions of the fre- 
quency u;, and 

^(O,^))--^^^^^^. (24) 



The asymptotically constant value of the optical-like 
plasmon frequency is then implicitly given by 

l-4Q%(0, W2 (0)) 5s A- 1 ^y(G)sin 2 Qg • 5 3 ) = 0, 

(25) 

whereas the imaginary part of the retarded polarization, 
close to the second plasmon branch [uo ~ ^2(0)], to lead- 
ing order in g, reads 



(19) Imn(q, W +) « -ngsA- 1 



MCia) -^F(CkA) 



1 



47V^(^ 2 (0)+4A-aA) 2 



x ^2 yrfhizij - yij)(zhk + yhk)qjQk5 - ^2(0)) (26) 

ijhk 

In particular, it follows that the spectral weight of Imll 
close to 002(0) decreases as ~ g 4 , as q — > 0, rather than as 
~ <? 3 / 2 , as is the case for the acoustic-like plasmon mode, 
Eq. ( [22] ) . This justifies the reduced spectral weight asso- 
ciated with the second plasmon branch at small wavevec- 
tor in Fig. [I] In the case of graphene on a dielectric sub- 
strate (e r > 1), inspection of Eqs. (24) and (25) yields a 
reduction of ^(O). 



III. EFFECT OF STRAIN ON THE PLASMON 
DISPERSION RELATION 

We now turn to consider the effect of strain on the 
plasmon dispersion relation. As in Refs. 25 28, applied 
uniaxial strain can be modeled by explicitly consider- 
ing the dependence on the strain tensor e of the tight- 
binding parameters tg = t(8g) through the vectors Si 
connecting two NN sites (£ = 1,2,3). A linear depen- 
dence of 81 on e is justified in the elastic limit. Such 
an assumption is however quite robust, due to the ex- 
treme rigidity of gra phene^, and is supported by ab ini- 
tio calculation J 19 | 21 l 

Below, the strain tensor e will be parametrized by a 
strain modulus e, and by the angle between the direc- 
tion of applied strain and the x axis in the lattice coordi- 
nate system. Specifically, one has = [resp., = 7r/6] 
for strain applied along the armchair [resp., zig-zag] di- 
rection. 

The main effect of applied strain on the band disper- 
sion relation is that of shifting the location of the Dirac 
points ±ko from their position ±K at e = 0. While va- 
lence and conduction bands vanish linearly as q —> ±ko 
for moderately low applied strain, such an approximation 
breaks down at a critical value of the strain modulus e, 
depending on the direction 6 of applied strain, when ±ko 
tends to either midpoint of the 1BZ border. This 
has been described in terms of an electronic topological 
transition (ETT), since it is accompanied by a change of 
topology of the Fermi line^. 

Fig. [2] shows the dispersion relation of the plasmon 
branches studied in Sec. |II C[ including LFE, along a 



FIG. 2: (Color online) Plasmon dispersion relation for doped 
graphene at finite temperature (fi = 1 eV, T = 3 K), including 
LFE, with strain applied along the = (armchair) direction. 
Strain increases (from left to right, from top to bottom) as 
£ = 0,0.075,0.175,0.275. 



FIG. 4: (Color online) Plasmon dispersion relation for doped 
graphene at finite temperature (ji = 1 eV, T = 3K), including 
LFE, with strain applied along the = tt/4 (generic) direc- 
tion. Strain increases (from left to right, from top to bottom) 
as e = 0, 0.075, 0.175, 0.275. 





FIG. 3: (Color online) Plasmon dispersion relation for doped 
graphene at finite temperature (fi = 1 eV, T — 3 K), including 
LFE, with strain applied along the = tt/6 (zig-zag) direc- 
tion. Strain increases (from left to right, from top to bottom) 
as e = 0, 0.075, 0.175, 0.275. 



FIG. 5: (Color online) Plasmon dispersion relation for doped 
graphene at finite temperature (/i = 1 eV, T = 3K), including 
LFE, with strain applied along the — (armchair) direction. 
Strain increases (from left to right, from top to bottom) as 
e = 0, 0.075, 0.175, 0.275. In each graph, different panels 
refer to cji(q) = wi(g,¥>q), with <p q = 0°, 20°, 40°, . . . 160°. 



symmetry contour of the 1BZ, for strain applied along 
the armchair direction (0 = 0), with increasing strain 
modulus (e = — 0.275). The low-frequency, 'acoustic' 
plasmon mode cJi(q) is not qualitatively affected by the 
applied strain. In particular, the dominant square-root 
behavior is independent with respect to the opening of a 
gap. On the other hand, one observes an increase of spec- 
tral weight associated with the high-frequency, 'optical' 
plasmon mode co>2(q) at small wavevectors. The overall 
flattening of the second plasmon branch over the sym- 
metry contour under consideration can be traced back to 
the strain-induced shrinking of both valence and conduc- 
tion bands. We also note the formation of a gap between 
^(q) and the pseudo-plasmon mode corresponding to a 
logarithmic singularity in Q°(q, oS) at q = K. 

A qualitatively similar analysis applies to the case 
of strain applied along the zig-zag direction (9 = 7r/6, 



Fig. [3J, and for strain applied along a generic direction 
(6 = 7r/4, Fig. [3]), with ^(q) dispersing more weakly as 
the strain increases. 

Finally, we turn to study the q-dependence of the low- 
frequency, 'acoustic' mode cji(q) = LJi(q, <Pq) under ap- 
plied strain, where q = |q| and <^ q denotes the angle be- 
tween q and the x axis. Fig. [5] shows then the dispersion 
relation of the lower plasmon branch as a function of q for 
several values of y? q , for increasing strain applied along 
the armchair direction (0 = 0). While the overall square- 
root shape uji ~ uJxJqa^ Eq. (21), is maintained in all 



cases, one observes a stiffening of such plasmonic mode 
with increasing strain and a maximum of the coefficient 
uji, Eq. (20), when <^ q — w 90°. The same descrip- 
tion qualitatively applies also to the cases of strain ap- 
plied along the armchair (0 = 7r/6), and along a generic 



7 



(6 = 7r/4) direction. Such a behavior can be jus tified 
analytically in the limit of no LFE (cf. Sec. II D), and 
corresponds to the strain dependence obtained for the 
optical conductivity^. Indeed, from Eq. (20), one may 



notice that all the strain dependence is contained in the 
modulus square of the quasiparticle dispersion relation of 
the conduction band at the Fermi level, |V q £?k2/&|- One 
finds 



tii oc |V q £ q2 | = 



cos 2 ((^ q -^) sin 2 ((/? q - 77) 



A 2 



£ 2 



1/2 



(27) 

where A and B denote the semiaxes of the constant en- 
ergy ellipse^, 



A- 2 
B~ 2 



with 



7 



3^ 



a\t{ + ti + ti), 



and 



cos(2?7) 
sin(2?7) 



a Va 2 + (3 2 



(28a) 
(28b) 

(29a) 

(29b) 
(29c) 

(30a) 
(30b) 



It follows that til attains its maximum values whenever 
(fq — rj = 7r/2 (modulo 7r), and its minimum values when- 



ever y? q — r\ = (modulo 7r). It turns out that rj = 6 
in the zig-zag and armchair cases (cf. Fig. [5|, whereas 
r] ~ in the generic case. 



IV. CONCLUSIONS 

By studying the electronic polarization, we have de- 
rived the dispersion relation of the plasmon modes in 
graphene. Besides including electron-electron correlation 
at the RPA level, we have also considered local field ef- 
fects, that are specific to the peculiar lattice structure un- 
der study. As a consequence of the two-band character of 
the electronic band structure of graphene, we find in gen- 
eral two plasmonic branches: (1) a low-energy, acoustic- 
like one with a square-root behavior at small wavevectors, 
and (2) a high-energy, optical-like mode, weakly dispers- 
ing at small wavevectors. This is generic to two-band 
systems, and might apply to other two-band systems as 
well, such as MgB2, and is analogous to the collective 
modes in two-dimensional quantum wells. We also find 
an intermediate energy pseudo-plasmon mode, associated 
with a logarithmic (viz. non power-law) divergence of the 
polarization, which can be related to an interband tran- 
sition between the Van Hove singularities in the valence 
and conduction bands of graphene, and can be identified 
with a 7r — > 7r* transition. We have next studied, both 
analytically and numerically, the dependence of the plas- 
mon branches on applied strain. While the square-root 
character of the low-energy mode at small wavevector is 
robust with respect to applied strain, we find a nonmono- 
tonic stiffening as a function of the wavevector direction, 
the maximum steepness occurring roughly when the lat- 
ter is orthogonal to the direction of applied strain. We 
have also studied the influence of applied strain on the 
high-energy, optical-like plasmon branch. 
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